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Abstract 



We determine non-perturbatively the anomalous dimensions of the sec- 



ond moment of non-singlet parton densities from a continuum extrapola- 
tion of results computed in quenched lattice simulations at different lattice 
spacings. We use a Schrodinger functional scheme for the definition of the 
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renormalization constant of the relevant twist-2 operator. In the region 
of renormalized couplings explored, we obtain a good description of our 
data in terms of a three-loop expression for the anomalous dimensions. 
The calculation can be used for exploring values of the coupling where a 
perturbative expansion of the anomalous dimensions is not valid a priori. 
Moreover, our results provide the non-perturbative renormalization con- 
stant that connects hadron matrix elements on the lattice, renormalized 
at a low scale, with the experimental results, renormalized at much higher 
energy scales. 
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1 Introduction 



The accurate knowledge of hadron parton densities is an essential ingredient for 
the experimental tests of QCD at accelerator energies. Their normalization is 
usually obtained from a fit to a set of reference experiments and used for pre- 
dicting the behaviour of hard hadron processes in different energy regimes. The 
energy evolution and the relative normalization of the hard processes are pre- 
dictable within renormalization-group-improved perturbation theory, while the 
calculation of the absolute normalization needs non-perturbative methods. These 
methods may provide not only a check of the non-perturbative aspects of QCD 
but also some information that could help in fixing the values of parton densities 
at values of the Bjorken x larger than 0.5 where the experimental information, 
especially for the gluon distribution, is scarce |TJ]. Besides, they could also help 
in establishing the domain of validity of a perturbative evolution that appears 
in some cases to be compatible down to scales where one would expect higher 
orders and power corrections to inelastic processes to take over. A study of these 
non-perturbative effects in moments of hadronic structure functions can be done 
within the lattice approach to QCD |2|, |3|, £3|. |J. 

In this paper we determine the non-perturbative anomalous dimensions of the 
second moment of non-singlet parton densities. The anomalous dimensions are 
extracted from a continuum extrapolation of results obtained in quenched lat- 
tice simulations at different lattice spacings. The definition of the renormaliza- 
tion constant and of the matrix element of the relevant twist-2 operator is done 
within the finite- volume Schrodinger functional (SF) scheme || [7], [8], |§. Non- 
perturbative renormalization constants within this scheme have been obtained 
by the ALPHA collaboration for local fermion bilinears, and in particular for the 
scale-dependent renormalization constant Zp entering the computation of the 
running quark mass [TO]. Here we follow ref. [IT], where the perturbative cal- 
culations needed to relate the SF scheme to more common schemes, such as the 
modified minimal subtraction scheme, have been carried out and some general 
considerations on the onset of lattice artefacts have been presented. 
In the region of couplings that we explore, we find a good agreement between 
our data and a three-loop approximation for the anomalous dimensions in the 
SF scheme. As a consistency check, we successfully compare the running of 
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our renormalization constant with the one expected from a fit to the anomalous 
dimensions. 

In section 2 we recall some basic facts about the Schrodinger functional scheme, 
we define the matrix element of the twist-2 operator and we outline the general 
strategy of the calculation. Section 3 contains the details of our numerical sim- 
ulations and the results for the continuum extrapolated values of the ratio of 
renormalization constants, i.e. the step scaling functions, calculated at renormal- 
ization scales differing by a factor 2. In section 4 we extract from these numerical 
results a fit to the renormalized coupling constant dependence of the anomalous 
dimensions of the operator. Section 5 is devoted to the conclusions. 

2 The Schrodinger functional scheme and the 
twist-2 operator 

The Schrodinger functional has been discussed extensively in the literature (see 



12[ |T3"1| for reviews); it represents the (quantum) amplitude for the time evolu- 
tion of a field configuration between two predetermined classical states at some 
time boundaries. It takes the form of a standard functional integral with fixed 
boundary conditions. It has been shown that its renormalizability properties are 
the same as those of the theory with infinite time extension, modulo the possible 
presence of a finite number of boundary counterterms J7|. In QCD with fermions, 
it can be written as: 



Z[C'J,p';C,p,p] = / V[U]Vty]V$]e 



-S[U,i/>,i/>] 



where C ,C and p',p',p,p are the boundary values of the gauge and fermion 
fields respectively. In the following discussion, the classical boundary gauge field 



will be set to zero. According to refs. || [TJ|], expectations values may involve 
the response ( to a variation of the classical Fermi field configurations on the 
boundaries: 



dp(x) 5p(x) 
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In the continuum, moments of non-singlet structure functions are related, through 
the operator product expansion, to hadronic matrix elements of local twist-2 
operators of the form: 



9 ^ n = {\T V^h^ B M2 • • • S M + trace terms , (4) 

where _D M is the covariant derivative, { "indices" } means symmetrization. 
The twist is defined as the difference between the engineering dimensions of the 
operator and its angular momentum. All listed operators belong to irreducible 
representations of the Lorentz group. 

On the lattice, the discretization of the covariant derivative can be done in a 
standard way: 



\J^{x) = -[Up(x)if>(x + ajl) - ip(x)] 

1 a (5) 

vjX^) = - U^x - ajl) 1 ip(x - aji)} . 

(X 

In this paper, we concentrate on the calculation of the second moment to which 
we associate the following irreducible operator that is multiplicatively renormal- 
izable: 



0\ 2 = ^7{l B 2} y V- ( 6 ) 

We define the SF matrix element of the second moment of the non-singlet parton 
densities by the observable: 



/ a (*o) = foM = -a 6 ^e^- z )(^(x)7[i D 2] ^x)C(y)r±r 3 C(z)> (7) 

where the contraction of the classical fields is non-vanishing if the matrix T sat- 
isfies: TP_( + ) = -P+(-), where P-{+) = |(1 ± 7o) and p is the momentum of the 
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classical field sitting on the boundary. The observable can be seen as the operator 
matrix elements between the vacuum and "p"-like classical states sitting at the 
T = boundary. 

We take the limit of massless quarks, which in the numerical simulations can be 
monitored via axial Ward identities. In the SF framework it is possible to work at 
zero physical quark mass because a natural infrared cutoff to the Dirac operator 
eigenmodes is provided by the time extent of the lattice. 

The matrix element of the operator for the second moment involves two directions: 
one of them is given by the contraction matrix T, i.e. from the polarization of 
the vector classical state: 



and the other one from the momentum p of the classical Fermi field at the bound- 
ary. For the tree level, this choice gives a non-vanishing matrix element in the 
massless quark limit, where we evaluate our correlations. The tree-level correla- 
tion can be easily calculated and reads: 



/f(*o) = ^ [(-<A) (M_(p)e-^ - M + (p)e- 2 ^--))] , (9) 



r = 7 2 , 




where 



P M = (2/a) sin(ap M /2), 



Pi* = (!/«) sin(a^) , 



(10) 



M(p) = m + -ap 2 , 



M ± = M(p) ± ip , 



(11) 




(12) 



R{p) = M{p) {1 



e -Mp)T|_^ o | 1+e -2a,(p)T| 



(13) 
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We have chosen as convention r = 72 and p = (pi,0, 0). Furthermore, we will 
always work on hypercubic lattices with T = L. 

We may define an unnormalized observable as the ratio of the correlation function 
/2OE0) divided by its tree- level expression 

with 7] < 1 and typically r\ — 1/2 or 77 = 1/4 (see below). The normalized 
observable is defined by removing the renormalization constant of the classical 
boundary sources (. Following refs. [16| , this is represented by the quantity 



Z(L) = , (14) 



called fi, also normalized by its tree-level expression. We thus arrive at our 
definition of the renormalization constant for the twist-2 non-singlet operator: 

Z{L) = ^=L . (15) 

We denote by L the physical lattice size and by a the lattice cut-off. The observ- 
ables on the lattice depend upon the ratio L/a, the value of the momentum api, 
and the point x^/a (or equivalently the value of 77) of the insertion of the operator 
in time. We note that, in fact, time translation symmetry is broken by the fixed 
boundary conditions. In the following, both p and Xq will be taken proportional 
to 1/L and L, respectively, so that the only external scale in the problem is L. 
In the expressions for Z, Z and fi the lattice spacing units are understood, and 
we will make them explicit only when needed. 

The operator needs renormalization to be finite in the continuum limit: we define 
the renormalization constant such that the operator matrix element is equal to 
its tree-level value at fi = 1/L. At one-loop, we obtain: 



R (n) = (1 - £ 2 Z«(l/a/i))O bare (a) 

(16) 

{l-g 2 ZW{l/a^){l + g 2 Z^{L/a))0^ 



where O tree is the result at zero coupling. Beyond perturbation theory, the renor- 
malized operator is defined by: 



R (fj) = Z(l/an)- 1 O baxe (a/L) (17) 
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with Z(L/a) defined by: 



O bare (a/L) = Z(L/a)O tiee , (18) 

where we have explicit the lattice cut-off dependence. The study of the scale 
dependence of the renormalization constant is then equivalent to that of the 
dependence upon the lattice size L, provided that, as anticipated, the external 
variables, upon which the matrix element may depend, scale like the basic length 
L in physical units: we hence take 

Pi = 2tt/L 

xq = r]L, (19) 

with 7] = 1/4 or 1/2. 

Following ref. fllfl , the renormalization constant up to one-loop, Z^(L/a), in the 
continuum can be parametrized by: 



Z« {L/ a) = B + C \n(L/a) , (20) 

with 

B = 0.2635(10) [q = 1/4] 

B = 0.2762(5) [rj = 1/2] . (21) 

The central goal of the calculation consists in obtaining the continuum limit of 
the so-called "step scaling function" <t z , describing the change of Z when the 
renormalization scale, proportional to the physical box size L, is varied by a 
factor of s, i.e. 

Z(sL) = a z (f(L))Z(L) , (22) 

at a fixed value of the running coupling g(L) renormalized at the scale L. The 
quantity o expressed in terms of the renormalized coupling constant does not 
depend anymore on the lattice cut-off and has a finite continuum limit. From the 
jump for a finite renormalization scale we can extract the one for an infinitesimal 
variation governed by the corresponding anomalous dimension. 



7 



On the lattice, we can only obtain the continuum quantities by extrapolating the 
lattice step scaling functions obtained, at fixed lattice spacing, from the ratio of 
renormalization constants computed on lattices with N and sN points. In order 
to extrapolate to the continuum, we perform several simulations at increasing 
values of the bare coupling constant and decreasing values of the lattice spacing 
by increasing the number of lattice points N, so that the physical volume L = 
Na remains constant. The latter condition is monitored by the value of the 
renormalized coupling constant in the SF scheme, where the renormalization scale 
is given by the physical lattice extent L. 

We have used the values of bare couplings and volumes corresponding to a fixed 
ot(L) determined by the ALPHA collaboration in their study of the running mass. 
In general, reaching large values of the renormalized coupling, and therefore large 
physical volumes, with a limited number of points, implies dealing with increasing 
lattice artefacts. 

We extrapolate our lattice step scaling functions to the continuum limit and 
compare the dependence of the continuum step scaling function on the coupling 
with the expression obtained from a perturbative expansion of the anomalous 
dimensions and of the /3-function in the renormalized coupling. 
More explicitly, the anomalous dimensions in a scale-dependent regularization 
scheme (momentum subtraction or SF scheme) are defined from the variation of 
the renormalization constant to an infinitesimal variation of the scale as: 



dlog(Z(/x)) 



dlog(/x) 



7(<7>)) (23) 



with g 2 (/i) satisfying: 



<hr(f') _ „ ,2 



• (24) 



dlog(jt) 

In our simulations we will always choose s = 2 in eq. fl22"|). The ratio Z(2L)/Z(L) 
can be expressed as: 
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Z(2L) r f 2L 

= exp^ 

^ (2L) , r 2/T ^dlog(Z(L')) dlog(L') 



- expj^ d[log(Z(L'))] 



For example, by inserting a perturbative expression to three-loop for gamma 

7 (s 2 (/i)) = 7 oS 2 (/i) + 7i^ 4 (/i) + 72S 6 (/i), (26) 

and for beta 

M^)) = A><7 4 M + (3i9 6 (v) + /WM, (27) 
we get for the finite scale jump: 

\og(Z(2L)/Z(L)) = F{g\2L)) - F(g\L)), (28) 

where 



F ^ = |Hh]!ot)-wH + 5' (,) (29) 

+ -§-(log(l + A/Ajx + A/ftx 2 ) - A/A>J(a;)) (30) 



with 



/(x) = (2/VA)arctg ^ l/A)a; ^ /A)X ), (31) 
A = 4/3 2 //3 - (A/A)) 2 , (32) 



for A > as in our case. 
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3 The step scaling functions 



In this section we present our numerical results for the non-perturbative eval- 
uation of the step scaling function. We used normal Wilson fermions without 
any improvement and worked on even-sized lattices ranging from 6 4 to 32 4 . We 



employed SSOR |20[ preconditioning and a BiCGstab solver |2T| for all necessary 
matrix inversions. The gauge fields were generated with a hybrid of heatbath 
and over-relaxation updates. We normally performed 12 to 16 over-relaxation 
steps per heatbath update and 20 to 50 iterations in between measurements for 
the 16 4 and the 32 4 lattices, respectively. All errors quoted below are computed 
using a jack-knife method. We explicitly checked by combining the jack-knife 
method with a binning procedure, that there exist no noticeable autocorrelation 
times for our observables. The statistics of our data are 300 measurements for the 
24 4 and the 32 4 lattices and reach 500 measurements for the smaller lattices. To 
complete the specification of the numerical simulations we performed, we finally 
give some Schrodinger functional specific parameters ||19|| : we set 9 = and the 
improvement coefficient Ct(go) to its one-loop value in order to cancel most of the 
extra O(a) corrections which would be absent in the pure gauge theory and are 
introduced through the Schrodinger functional boundary conditions. In addition 
we used trivial background gauge fields. 

We start the discussion of our numerical results by showing in fig. [I] a typical 
correlation function, /^(^o) from eq.(^), computed on a 16 4 lattice. The signal 
can be followed up to large distances and only for, say, Xq > 3L/4, the correlation 
function becomes too noisy. A similar qualitative behaviour of our correlation 
functions was found for other lattices, too. For the definition of Z, we can use 
Xq = L/4 or Xq = L/2, where the errors of the correlation functions are reasonably 
small. However, the continuum extrapolation of the results at Xq = L/2 appears 
more affected by lattice artefacts, as already observed in ref . []TTJ ; therefore, in the 
following, we will only present results for the case xq = L/4, i.e. we will choose 

V = \ (33) 

from now on. We remark that the case L/a = 6 for xq = L/4 is obtained by 
an interpolation from the points at xq = (L ± 2a)/ 4 that cancels leading lattice 
artefacts. 
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Figure 1: An example for a correlation function /^(^o) m eq. (0), for a 16 4 
lattice. Parameters are /3 = 7.0197 and k = 0.144987. 



11 



In order to obtain the running of the Z(/i), eq. (lq), we computed the step scaling 
functions for Z(L): 

Z(2L) . , 

" = W' (34) 

for /i 



and for the renormalization constant itself, 

^ = w (36) 

We first discuss er^ and erg. We have calculated cr^ and erg for five values of 
the coupling in the Schrodinger functional scheme renormalized at scale L. All 
results are reported in table 1. There we also give the values of (3 taken for 
the continuum extrapolation of the above-mentioned ratios. We recall that the 
values of f3 are chosen such that when varying the number of lattice points N, 
the physical length L of our box stays constant. 

In fig. H] we show our data for erg. The different figures are labeled by the value 
of the running coupling constant g 2 (L), taken at scale L, in the Schrodinger 
functional scheme. For erg we see a marked curvature of the data when plotted as 
a function of a/L. We tried to fit the data using a linear and quadratic (in a/L) 
ansatz. For the linear fit, we used the three data points with smallest values of 
a/L, whereas for the quadratic fit we used all data points. For most cases, the 
value of x 2 P er degree of freedom (d.o.f.) was much better for the fit using a 
quadratic form. It seems that the use of non-improved Wilson fermions induces 
large lattice artefacts that appear linear and quadratic in a/L. We think that 
these large lattice artefacts are due to the non-zero momentum that we are using 



for computing our observable. This is in accord with the remarks in [IT]. As can 
be seen in table 2, the values of erg, extrapolated to the continuum, suffer from 
the large error induced by the use of the quadratic extrapolation. 
The situation is somewhat different for the continuum extrapolation of cr^, shown 
in fig. (3). Here we find that a linear fit in a/L, again leaving out the value of 
er/j for the scaling from L/a = 6 up to L/a = 12, always gives a better or 
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7.3686 


0.0625 


0.8981(21) 


0.9095(44) 


7.1190 


0.0833 


0.9041(18) 


0.9285(45) 


6.7860 


0.1250 


0.9020(28) 


0.9541(63) 


6.5512 


0.1667 


0.8824(27) 


1.0067(77) 


7.0203 


0.0625 


0.8762(24) 


0.9012(55) 


6.7750 


0.0833 


0.8867(24) 


0.9172(55) 


6.4527 


0.1250 


0.8792(21) 


0.9622(55) 


6.2204 


0.1667 


0.8729(22) 


1.0451(68) 



Table 1: The lattice step scaling functions used for the continuum extrapolations 
given in table 2 and table 3. 
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Figure 2: Continuum extrapolation of oz to obtain the step scaling function for 
all values of the renormalized coupling considered. The solid line is a quadratic 
fit to the data. The values of g 2 , labelling the graphs, are taken at scale L. 
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compatible y 2 /d.o.f. as compared to a quadratic fit taking all data points. Since 
/i is computed at zero momentum, the different behaviour of <jf 1 as compared to 
Oz can be traced back again to the lattice artefacts associated to the non-zero 
momentum present only in the unnormalized constant Z . 

From the above discussion it follows that and 0% have to be extrapolated 
differently to the continuum limit, and we decided to perform this extrapolation 
independently for the step scaling functions. Such a strategy is certainly justified, 
given the fact that we did not find any autocorrelation time in our data for fi and 
Z, and that there was only little correlations between the two quantities. (Re- 
member that f\ is computed at zero momentum, while Z needs a non-vanishing 
momentum.) We then finally compute the step scaling function for the renormal- 
ization constant of the twist-2 non-singlet operator as the ratio of the individual 
continuum extrapolations of a 2 and a^. 

We summarize our results for all values of the different step scaling functions 
in table 2, together with the corresponding values of g 2 (L) in the Schrodinger 
functional scheme. We summarize the continuum values of az in table 3. 





°h 


x 2 


&z 


x 2 


Type of fit 


1.8811 


0.919(19) 


3.95 


0.917(8) 


0.19 


quadratic 


2.1000 


0.894(23) 


1.26 


0.891(9) 


0.20 


quadratic 


2.4484 


0.931(23) 


4.24 


0.891(9) 


1.01 


quadratic 


2.7700 


0.888(29) 


1.66 


0.854(13) 


0.01 


quadratic 


3.48 


0.913(30) 


0.31 


0.855(12) 


6.44 


quadratic 


1.8811 


0.912(7) 


2.61 


0.931(3) 


1.53 


linear 


2.1000 


0.897(9) 


1.17 


0.920(4) 


3.60 


linear 


2.4484 


0.886(9) 


6.90 


0.913(4) 


0.03 


linear 


2.7700 


0.866(10) 


0.52 


0.897(5) 


3.64 


linear 


3.48 


0.837(11) 


0.39 


0.880(5) 


10.2 


linear 



Table 2: The values of the step scaling functions are given in the continuum. The 
running coupling is computed in the SF scheme. 
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Figure 3: The continuum extrapolation of Of x . The solid line is a linear fit to 
the data. The values of g 2 , labelling the graphs, are taken at scale L. The error 
bars in this case are smaller than the symbol sizes. 
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cz 


1.8811 


1.006(12) 


2.1000 


0.993(15) 


2.4484 


1.006(16) 


2.7700 


0.985(20) 


3.48 


1.021(21) 



Table 3: The continuum values of the step scaling function for the renormalization 
constant of the twist- 2 non-singlet parton density. 



4 The non-perturbative anomalous dimensions 

The numerical results discussed in the previous section allow us to extract the 
non-perturbative values of the anomalous dimensions. We want to emphasize 
that, from now on, we are discussing results that are already extrapolated to 
the continuum in the, maybe, somewhat unusual Schrodinger functional scheme 
(instead of, say, the MS scheme). 

The extraction of the anomalous dimension through the step scaling function, in 
general requires the knowledge of the beta function to the appropriate accuracy. 
For the running coupling g 2 in the SF scheme, the beta function is well expressed 
by a three-loop formula up to g 2 = 3.5 |]19|. However, the range of renormalized 



couplings explored in our work (see table 2 and in particular the last two values 
g 2 (L) = 2.77 and g 2 (L) = 3.48) lead to values for g 2 (2L) outside the domain 
of validity of the three-loop parametrization. Therefore, the extraction from 
eq. fl2"8]), valid to three-loop only, of the two- and three- loop coefficients of the 
perturbative expansion of the anomalous dimensions in g 2 {L), can only be done 
safely for the first three values of g 2 (L) simulated. 

To describe our data we have performed fits to the step scaling functions using 
eq. fl28|), i.e. by expanding ■y(g 2 ) and f3(g 2 ) up to three-loop order, eqs. (p6j), (p7|). 
For the fit for az we have fixed the one- and two-loop contribution by setting 70 
and 71 to their perturbative values. The latter are extracted from a conversion 



from the modified minimal subtraction scheme results JT8| to the SF scheme. It 
takes into account both the different operator renormalization constant B in 
eq. (m 



17 



B = 0.2635 (SF) 

B = 0.0108 (MS) (37) 



and the relation between and p| F |T 



o 9 1 .2556 ^ /r* \ 

%s = ^sf + -J^y9sF- (38) 

The perturbative two-loop coefficient to be expected is: 

7l SF = -0.0268 [fi- 1 = L] , (39) 

where = L indicates the scale at which the SF coupling is renormalized, to 
be compared with 

7 p = 0.0039. (40) 

The fit to our data of <Jz is shown in fig. |j. We find for the three-loop coefficient 
72: 

7 f = 0.0034(7) [/i- 1 = L] . (41) 

We also plot in fig. f| the analytical form of the step scaling function when we 
truncate 7(g 2 ) to one- and two- loop only, while always keeping the three- loop 
expansion of I3{g 2 ). We remind that the fit is based on the first three points only, 
although we show in fig. |] also the remaining data points. 

The corrections due to higher loops are large and oscillating in sign, which might 
be a signal of an unfortunate choice of the expansion parameter renormalized 
at the scale /i -1 = L. Therefore, we have repeated the fit for the step scaling 
function by using a different scale for the running coupling, g^ F (L/4) instead 
of p| F (L). This opens the possibility of including in the fit the largest ggp^) 
points that correspond in the case of p| F (L/4) to still moderate values where the 
three-loop parametrization is valid. The choice of the scale was motivated by the 
identification with the value of Xo = L/4 where the operator is inserted and we 
did not attempt to optimize such a scale. 
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Figure 4: The dependence of o z on the running coupling g%${V) in the 
Schrodinger functional scheme taken as r\ — 1, i.e. the scale to be L. 
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We show a fit to our data, employing the scale L/A in the upper graph of fig. |^. 
As above, the 70 and 71 were fixed to their perturbative values. The new value 
for 7l F obtained from eq.(39) after the change of scale is: 

7l SF = -0.0181 [fi- 1 = L/A] , (42) 

and we obtain from our fit for the three- loop coefficient: 

7 f = -0.005(3) [fi- 1 = L/A] , (43) 

with a x 2 P er degree of freedom of 0.7. Figure |5] shows that, with the choice of 
a smaller scale and of correspondingly smaller values for the running couplings, 
the relevance of higher-loop terms decreases and a three-loop expression appears 
to be safe. 

We have also made a fit to the anomalous dimensions of the wave function renor- 
malization constant v/i : i n this case we can only fix 70 to its perturbative value 
70 = 0.05066. Fitting both the two- and three-loop coefficients leads to large 
errors in the fitted values. We therefore attempted only a fit of the two-loop co- 
efficient which gave a reasonable value of x 2 — 0.1. Using the coupling g^ F (L/A), 
the fit turns out to be very stable when three, four or five points are included in 
the fit. The value for the two-loop coefficient obtained is 

7 fK = 0.030(3) , (44) 

The result for p| F (L/4) is plotted in fig. (lower graph). In general, the size- 
able two- and three-loop contributions for Z and f\ are mainly due to the large 
constants Bq appearing in the definition of the renormalization constants. When 
such a constant is small, as in the case of the unnormalized Z, the two-loop co- 
efficient is rather small, of the order of the MS coefficient. We attempted a fit to 
a 2 and find: 

7 f = 0.003(3) , (45) 
indicating a dominant one-loop running. 

5 Conclusions 

We have studied the non-perturbative running of the average momentum of non- 
singlet parton densities in the SF scheme in the region of renormalized o>sf ranging 
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g§p(L/4) 



Figure 5: The dependence of oz and cr^ on the running coupling in the 
Schrodinger functional scheme evaluated at scale L/4, ^| F (L/4). 
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Figure 6: The dependence of Z on the scale /i as compared to the three-loop 
running, where /xo is the scale at which Z is normalized to unity. 

from 0.1 to 0.2. From the step scaling functions we have extracted a three- 
loop parametrization of the anomalous dimensions. In turn, from the knowledge 
of the three-loop anomalous dimensions, we can calculate the running of the 
renormalization constant Z . As a check, this running is compared in fig. ^| with 
our results for the step scaling functions. The running in the SF scheme appears 
to be rather slow for the energy ranges that we have explored. The errors of the 
data obtained in our numerical simulations are sizeable. 

The values of a s covered in this work, according to ref. ||19|| , correspond to the 
energy range from 1 to 10 GeV. Higher values of the renormalized couplings can 
be explored, and a higher accuracy for the continuum extrapolation reached if 
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leading lattice artefacts, mainly related to the lattice momentum quantization, 
are removed, for example by using the non-perturbatively improved clover ac- 
tion and, possibly, improved operators. This calculation allows us to connect 
the continuum limit of a lattice evaluation of the hadron matrix element of the 
average momentum of non-singlet parton densities, renormalized at a fixed low- 
energy scale, with high energy experimental results, without an a priori unknown 
systematic error deriving from the truncation of the perturbative series for the 
anomalous dimensions. 

ACKNOWLEDGEMENTS. We thank H. Wittig for sending us,prior to publica- 
tion, results for the values of /3 used in this work to perform the matching. We 
are most grateful to M. Liischer, R. Sommer and P. Weisz for critical discussions. 

References 

[1] J. Huston et al., Phys.Rev. D58 (1998) 114034. 

[2] G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas, Nucl. 
Phys. B445 (1995) 81. 

[3] M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz 
and A. Schiller, Nucl. Phys. B472 (1996) 309. 

[4] C. Best, M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, A. 
Schafer, G. Schierholz, A. Schiller and S. Schramm, Phys. Rev. D56 (1997) 
2743. 

[5] M. Gockeler, R. Horsley, H. Oelrich, H. Perlt, D. Petters, P.E.L. Rakow, 
A. Schafer, G. Schierholz and A. Schiller, preprint DESY-98097, hep- 
lat/9807044. 

[6] M. Liischer, R. Narayanan, P. Weisz and U. Wolff, Nucl. Phys. B384 (1992) 
168. 

[7] S. Sint, Nucl. Phys. B421 (1994) 135. 

[8] M. Liischer, S. Sint, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B491 
(1997) 232. 



23 



[9] K. Jansen, C. Liu, M. Liischer, H. Simma, S. Sint, R. Sommer, P. Weisz and 
U. Wolff, Phys. Lett. B372 (1996) 275. 

[10] S. Capitani, M. Guagnelli, M. Liischer, S. Sint, R. Sommer, P. Weisz and 
H. Wittig, Nucl. Phys. B (Proc.Suppl.) 63 (1998) 153. 

[11] A. Bucarelli, F. Palombi, R. Petronzio and A. Shindler, CERN preprint, 
CERN/TH-98-251, hep-lat/98-08-005. 

[12] R. Sommer, DESY preprint, DESY 97-207, hep-lat/9711243. 

[13] M. Liischer, DESY preprint, DESY 98-017, hep-lat/9802029. 

[14] M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B478 (1996) 365. 

[15] M. Liischer, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B491 (1997) 344. 

[16] S. Sint and P. Weisz, preprint, MPI-PhT/98-50, hep-lat/9808013. 

[17] S. Capitani and G. Rossi, Nucl. Phys. B433 (1995) 351. 

[18] E.G. Floratos, D.A. Ross and C.T. Sachrajda, Nucl. Phys. B129 (1977) 66. 
Erratum: B139 (1978) 545. 

[19] M. Liischer, R. Sommer, P. Weisz and U. Wolff, Nucl. Phys. B413 (1994) 
481. 

[20] S. Fischer et. al, Comp. Phys. Comm. 98 (1996) 20; N. Eicker et. al., Nucl. 
Phys. B (Proc. Suppl.) 63 (1998) 955. 

[21] A. Frommer, V. Hannemann, T. Lippert, B. Noeckel and K. Schilling, Int. 
J. Mod. Phys. C5 (1994) 1073. 



24 



